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Section I 

INTRODUCTION 


The Solid Rocket Booster (SRB) of the Space Shuttle is expected to descend 
and impact into the ocean shortly after launch. Current plans call for the 
recovery of the SRB^s from the ocean. In order to estimate the probable loss 
or damage to the SRB’s in the recovery process, the following information is 
needed : 

• Ocean current and direction distribution 

• One-kilometer altitude winds 

• Wave height distribution 

• Wave slope distribution 

• Vertical wave velocity distribution 

• Horizontal wave velocity distribution. 

This information is required at the two locations presently being considered 
for ocean entry and recovery operations. These locations are the Cape Kennedy, 
Florida, Atlantic coastal waters bounded by 24 to 30 degrees north latitude, 

75 to 80 degrees west longitude; and the Vandenberg AFB, California, coastal 
waters, bounded by 31 to 33 degrees north latitude, 120 to 122 degrees west 
longitude. 

To improve our estimate of losses of and damage to the SRB, a simulation 
of the behavior of the ocean waves would be valuable. Such a simulation is 
important in order to assess the interaction of the SRB with the ocean and to 
calculate the stress loadings caused by the waves. In this study, only methods 
for simulating long-crested waves are developed. The simulation of more 
general cases remains for future work. 

Extensive work has been done in the study of ocean waves (References 7 
through 11), but substantially less work has been done in ocean simulations. 

No previous work, as far as we know, has been carried out using the control 
system model approach to ocean wave simulation. 
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This report also discusses and presents in tabular form the statistical 
distribution and sampling procedures to be used in simulation studies for the 
ocean currents (Section II) and for the 1-kilometer winds (Section III) . In 
Section IV the long-crested wave model is developed and discussed, the Pierson- 
Moskowltz spectrum is given, and the equations are normalized. In Section V, 
two simulation procedures are developed and described. The statistical distri- 
butions of the wave height, wave slope, and surface velocity components, as 
well as velocities below the surface, are given in Section VI. 
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Section II 


OCEAN CURRENTS 


2.1 GENERAL 

To simulate an ocean current for a study of the Solid Rocket Booster (SRB) 
splashdown and recovery, the joint probability distribution ^ of the 
current speed w and direction 0 are required. The term f^ ^ is given by the 
product of the marginal distribution of w, f^, and the conditional distribution 
of 0 given w, f^|^. 

The simulation of ocean currents requires the sampling of values of w and 

0 from f _ I f . The sampling procedure first samples w from the marginal distri- 
0 )w w 

bution, f ; then this value is used in the conditional distribution f ^ i to 
’ w’ 0 |w 

sample the value of 0 . 

Computer subroutines are available to generate uniformly distributed 

pseudo-random numbers, R^, between 0 and 1. A common method of random sampling 

from a distribution f (w) is to first randomly obtain R , then use a transfor- 

w 

mat ion to obtain the sampled value of w. This relationship is given by (see 
Reference 1) 


R 

w 


w 


fW 

f(w’) dw’ 
o 


( 2 . 1 ) 


where F 
w 


is the cumulative distribution of the ocean current speed. 


Thus, the procedure is to have the computer generate a value of R^, then 
solve Equation 2.1 for w. This is accomplished by using an interpolation scheme 
as follows. At discrete values of F^, [F^ = 0, .1, .2, ..., 1.0], the corres- 
ponding values of w are obtained. These values are inputted into a spline 
curve-fitting routine along with the values of the Initial and final derivatives. 


The routine fits cubic curves between adjacent input values, and it requires 
continuous values of the function and its first derivative at each point. 

Then for each sampled value of which is set equal to F(w), the appropriate 
sampled value of w is obtained by interpolation. 

The spline-fitting using the input data is available on the UNIVAC 1108 
as a subroutine SPLNl (Reference 2). To obtain interpolated values of w for 
given values of R^, use is made of the subroutine SPLN2 . This procedure will 
be illustrated in the following subsections. 

2.2 OCEAN CURRENT SPEED 

The cumulative distribution of the ocean current speed F needed for the 

^ w 

random sampling is given in Table 2-1 along with the required derivatives 
dw 

-TTr at F =0 and at F =1. These results were obtained from Reference 3. 
dF^ w w 

The general terms, "summer” and "winter" in the table refer to the months. 
May through October and November through April, respectively. The three 
oceanic areas are areas 22 and 26, which are off the Atlantic coast near 
Kennedy Space Center (KSC) , and an area which is off the Pacific coast near 
Vandenberg Air Force Base (VAFB). These areas are shown on Figure 2-1. 

In the simulation procedure used to sample randomly from the distribution 
f^, the data F^ and the corresponding values of w given in Table 2-1 are first 
inputted into the subprogram, SPLNl. Then a random number R^ is generated. 
Setting R^ equal to F^ and using the subprogram SPLN2, the interpolated value 
of w is obtained for use in the simulation. 

On Figure 2-2, the values entered into the spline routine are shown, as 
well as the interpolated values found by the spline routine for an illustrative 
case. 


2.3 CONDITIONAL OCEAN CURRENT DIRECTION 

Once the ocean current speed w has been sampled, the ocean current 
direction 0 is then randomly sampled for each given value of w. Similar to 
Equation (2.1), Equation (2.2) can be written as 
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OCEAN CURRENT SPEED AND DIRECTION DATA WERE PROCESSED FOR THREE DESIGNATED AREAS, NAMELY 
AN AREA ALONG THE LOWER CALIFORNIA COAST AND AREAS 22 AND 26 OFF THE FLORIDA COAST. THE 
CALIFORNIA AREA VARIES SEASONALLY FROM THE BOUNDARIES SHOWN HERE, WHEREAS AREAS 22 AND 26 
ARE UNCHANGING. 



LOCATION OF AREA ALONG THE LOWER 
CALIFORNIA COAST IN THE VICINITY OF 
VANDENBERG AFB (VAFB) 



LOCATION OF AREAS 22 AND 26 OFF THE 
FLORIDA COAST IN THE VICINITY OF 
KENNEDY SPACE CENTER (KSC) 


1 

I 

8 

! 


I 


Figure 2-1. LOCATION OF SRB IMPACT ZONES 






( 2 . 2 ) 



where is the conditional, cumulative distribution of 0 for a given value 

of w. Thus, the same procedure is followed for 0 as for w; namely, another 
random number, R , is generated and set equal to F i to obtain the value of 0, 

0 0 jW 

using the spline curve interpolation of the appropriate cumulative distribution. 


The conditional, cumulative distribution for 0 and the derivatives at the 
extrema of are given in Table 2-2 for use in the spline interpolation pro- 
gram. These are given for the same geographic regions and time periods as in 
subsection 2.2. "Summer** and **winter*' have the same measurable ranges as in 
Table 2-1. 


In Figure 2-3, an illustrative set of data points to be inputted into the 
spline routine is shown along with the interpolated values obtained from the 
spline routine. 


By generating two random numbers, R and R , a selection can be made from 

w 0 

the appropriate statistical distribution sample values of ocean current speed 
w and ocean current direction 0 for use in SRB ocean studies. A large number 
of sample values of w and 0 will satisfy the appropriate distribution of 


f 


0,w* 
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Table 2-2. CONDITIONAL CUMULATIVE DISTRIBUTION 
GIVEN OCEAN CURRENT SPEED (KNOTS), 


CUMULATIVE FREQUENCY 


LOCATION 

a* 

0 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

KSC 










AREA 22: 
SUMMER 
0<w<l .0 

0 

0 

0 

12 

45 

135 

164 

184 

197 

1.0<w 

0 

0 

0 

0 

16 

39 

151 

176 

186 

WINTER 










0<w<1.0 

0 

0 

0 

24 

86 

148 

173 

197 

202 

1 . 0<w 

0 

0 

0 

11 

123 

165 

184 

202 

206 

AREA 26: 
SUMMER 
0<w<l .0 

0 

0 

0 

0 

0 

0 

6 

19 

31 

l7o<w 

0 

0 

0 

0 

0 

0 

0 

3 

12 

WINTER 
0<w<l .0 

0 

0 

0 

0 

0 

0 

0 

0 

15 

1 .0<w 

0 

0 

0 

0 

0 

0 

0 

5 

16 

VAFB 










AREA 23: 

SUMMER 

0<w<0.9 

0 

0 

0 

63 

100 

125 

146 

164 

190 

0.9<w 

0 

0 

0 

62 

84 

103 

122 

141 

164 

WINTER 

0<v;<0.9 

0 

0 

19 

80 

103 

118 

133 

148 

165 

0.9<w 

0 

40 

75 

92 

103 

111 

119 

127 

136 

* a - ^ 

. 1^9 

|w " 

0 







** b = ^ 
^1, 


|w " 

X 0 

1 

.01 







OCEAN CURRENT DIRECTION (DEGREES), 
SPLINE CURVE INTERPOLATION 


.80 .90 1.00 b** 


213 251 360 9.0 

192 202 360 5.0 

217 253 360 9.0 

218 251 360 9.0 

97 140 360 6.0 

22 37 360 14.0 


32 72 360 11.0 

28 45 360 15.0 


232 272 360 8.0 

273 285 360 6.0 


190 246 360 10.0 

143 154 360 6.0 







Section III 


ONE-KILOMETER WINDS 


3.1 GENERAL 

For the purpose of simulating ocean conditions, knowledge of the joint 
distribution of wind direction and wind velocity at the 1-kllometer altitude 
is also required. As in subsection 2.1, the joint distribution of the wind 
speed V and the wind direction ((> are written as a product of the marginal 
distribution of v, f^, times the conditional distribution of ^ for a given 
As before 

(|),V (j) |v V 

3.2 ONE-KILOMETER WIND SPEED 

The cumulative distribution of v, F^, is represented in tabular and graphi- 
cal form in Table 3-1 and Figure 3-1. These values were obtained from Reference 

4. Table 3-1 gives the marginal, cumulative distribution values, F for the 

( dv 1 

^ j at F^ = 0 

ana at = x, ror use xn tne spxxne xnterpoxatxon. xo sampxe v from f^ we 
generate from a uniform distribution, set it equal to F^, then find v by 
interpolation . 


3.3 ONE-KILOMETER WIND DIRECTION 

Analogous to subsection 2.3, once the wind velocity has been sampled, the 
wind direction is randomly sampled. This can be done through the conditional, 
cumulative distribution of wind direction (j) for a given value of wind 

velocity v. The values needed for implementing the spline fit are given in 
Table 3-2 for the KSC geographical area and in Table 3-3 for the VAFB geographi- 
cal area. Examples of the cumulative distribution of data points used in the 
spline fit and the interpolated results are presented on Figure 3-2. To sample 

the wind direction <f) we generate a random number R from a uniform distribution, 

<P 

set it equal to where the v has been obtained previously, and interpolate 

the random value (f). 


The discussion in subsection 2.3 on sampling ocean currents is applicable 
here to the sampling of winds. 
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Table 3-1. MARGINAL CUMULATIVE DISTRIBUTION OF SCALAR WIND SPEED 
AT 1 KM (m/sec) 


LOCATION 

a* 

0 

.10 

.20 

CUMULATIVE FREQUENCY (Fy) 

.30 .40 .50 .60 .70 .80 

.90 

1.00 

b** 

KSC 

JAN 

0.4 

0 

2.7 

4.3 

5.6 

6.8 

7.4 

9.1 

10.4 

12.0 

15.0 

26.9 

5.7 

FEB 

0.3 

0 

2.7 

4.2 

5.5 

6.5 

7.5 

9.0 

10.9 

13.1 

16.3 

28.0 

3.5 

MAR 

0.0 

0 

1.5 

2.2 

2.8 

3.3 

3.6 

4.5 

5.5 

6.7 

9.5 

22.7 

3.0 

APR 

0.2 

0 

2.2 

3.6 

4.7 

5.7 

6.6 

7.5 

8.6 

9.8 

13.0 

24.0 

3.4 

MAY 

0.1 

0 

1.7 

2.6 

3.5 

4.3 

5.0 

5.9 

7.0 

8.1 

9.8 

18.6 

4.0 

JUN 

0.1 

0 

1.6 

2.4 

3.1 

3.8 

4.5 

5.2 

6.1 

7.1 

9.3 

25.0 

12.0 

JUL 

0.1 

0 

1.6 

2.4 

3.1 

3.8 

4.5 

5.2 

6.1 

7.1 

9.2 

24.9 

12.0 

AUG 

0.0 

0 

1.3 

2.2 

3.0 

3.7 

4.3 

5.0 

5.8 

6.7 

8.2 

20.9 

9.0 

SEP 

0.1 

0 

1.5 

2.5 

3.5 

4.4 

5.4 

6.5 

7.6 

9.0 

11.7 

32.5 

11.0 

OCT 

0.1 

0 

2.1 

3.5 

4.4 

5.2 

6.1 

7.1 

8.3 

9.7 

11.8 

19.9 

3.5 

NOV 

0.3 

0 

2.6 

3.8 

4.8 

5.6 

6.5 

7.5 

8.6 

10.0 

12.2 

23.5 

4.0 

DEC 

0.4 

0 

2.7 

4.0 

5.2 

6.2 

7.1 

8.2 

9.4 

10.7 

12.9 

22.9 

4.0 

VAFB 

JAN 

0.1 

0 

1.4 

2.6 

3.7 

4.7 

5.6 

7.1 

8.5 

10.2 

12.7 

23.7 

5.0 

FEB 

0.1 

0 

1.8 

3.1 

4.4 

5.6 

6.5 

7.7 

9.0 

10.7 

13.1 

20.4 

8.0 

MAR 

0.2 

0 

2.0 

3.1 

4.3 

5.4 

6.5 

7.7 

8.9 

10.1 

12.2 

18.7 

8.0 

APR 

0.3 

0 

2.5 

3.7 

4.8 

5.7 

6.5 

7.4 

8.4 

9.5 

11.2 

21.4 

12.0 

MAY 

0.1 

0 

1.5 

2.6 

3.5 

4.3 

5.0 

6.1 

7.3 

8.6 

10.3 

17.3 

8.0 

JUN 

0.1 

0 

1.5 

2.4 

3.4 

4.0 

4.6 

5.6 

6.7 

7.8 

9.5 

19.4 

12.0 

JUL 

0.1 

0 

1.1 

1.9 

2.5 

3.0 

3.5 

4.2 

4.9 

5.8 

7.3 

14.3 

8.0 

AUG 

0.1 

0 

1.2 

1.7 

2.0 

2.3 

2.6 

3.7 

4.9 

6.2 

8.2 

13.7 

6.0 

SEP 

0.1 

0 

1.0 

1.7 

2.1 

2.4 

2.7 

3.7 

4.8 

6.3 

8.2 

13.4 

6.0 

OCT 

0.1 

0 

1.0 

1.8 

2.8 

3.6 

4.3 

5.3 

6.3 

7.7 

9.8 

22.3 

14.0 

NOV 

0.1 

0 

1.5 

2.5 

3.5 

4.3 

5.2 

6.6 

8.0 

9.7 

11.9 

24.4 

13.0 

DEC 

0.2 

0 

1.8 

3.2 

4.4 

5.4 

6.5 

7.8 

9.2 

10.8 

13.2 

20.3 

8.0 


a = 


dv 


dF 


v'F=0 


** b = 


dv 


dF. 


F=1 


3-2 



WIND VELOCITY, 
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Table 3-2. CONDITIONAL CUMULATIVE DISTRIBUTION OF WIND DIRECTION (DEGREES), 
GIVEN WIND SPEED AT 1 KM (m/sec), AT KSC GEOGRAPHICAL AREA. 

FOR SPLINE CURVE INTERPOLATION 


KSC 

m/sec 

a* 

0 

.10 

.20 

.30 

CUMULATIVE FREQUENCY 
.40 .50 .60 .70 

.80 

.90 

1.00 

b** 

JAN 

0<v^5 

2 

0 

49 

113 

147 

166 

194 

225 

252 

279 

311 

360 

7 


5<v^l0 

6 

0 

61 

92 

127 

157 

195 

240 

264 

290 

311 

360 

8 


10<v<15 

9 

0 

79 

171 

197 

218 

232 

250 

267 

289 

311 

360 

9 


15<v^20 

34 

0 

189 

210 

236 

245 

254 

263 

275 

286 

300 

360 

5 


20<v^25 

31 

0 

173 

188 

197 

205 

210 

222 

237 

256 

275 

360 

12 


25<v 

29 

0 

170 

185 

190 

195 

203 

209 

236 

255 

275 

360 

12 

FEB 

0<v<5 

4 

0 

52 

97 

136 

161 

187 

228 

259 

279 

308 

360 

7 


5<v<10 

8 

0 

82 

130 

162 

182 

205 

230 

252 

275 

299 

360 

9 


10<v<15 

8 

0 

117 

172 

192 

211 

235 

259 

273 

285 

306 

360 

8 


15<v<20 

20 

0 

170 

192 

213 

238 

250 

261 

273 

285 

297 

360 

11 


20<v<25 

38 

0 

196 

208 

216 

228 

238 

245 

257 

269 

286 

360 

13 


25<v 

40 

0 

208 

213 

219 

224 

230 

236 

245 

258 

280 

360 

13 

MAR 

0<v<5 

2 

0 

21 

96 

127 

161 

180 

212 

239 

265 

307 

360 

7 


5<v^lO 

6 

0 

59 

115 

145 

175 

204 

235 

261 

279 

309 

360 

7 


10<v<15 

10 

0 

148 

188 

206 

219 

236 

255 

268 

286 

307 

360 

7 


15<v<20 

15 

0 

144 

193 

211 

224 

237 

251 

262 

278 

297 

360 

7 


20<v 

20 

0 

139 

199 

217 

230 

239 

248 

256 

270 

287 

360 

7 

APR 

0<v<5 

2 

0 

49 

87 

120 

151 

179 

204 

237 

274 

307 

360 

7 


5<v^l0 

6 

0 

64 

104 

126 

144 

175 

214 

244 

270 

285 

360 

7 


10<v<15 

11 

0 

81 

101 

121 

146 

207 

240 

264 

285 

312 

360 

7 


15<v<20 

0 

0 

175 

207 

220 

230 

238 

251 

262 

270 

285 

360 

7 


20<v 

0 

0 

210 

220 

230 

240 

252 

261 

267 

275 

292 

360 

7 

MAY 

0<v<5 

1 

0 

29 

65 

89 

117 

149 

178 

208 

248 

294 

360 

7 


5<v<10 

: 6 

0 

53 

82 

98 

121 

149 

178 

213 

242 

285 

360 

7 


10<vil5 

3 

0 

30 

47 

63 

87 

121 

196 

230 

254 

268 

360 

16 


15<v 

1 

0 

14 

34 

56 

221 

229 

236 

246 

259 

281 

360 

12 


* 




** 
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Table 3-2. CONDITIONAL CUMULATIVE DISTRIBUTION OF WIND DIRECTION (DEGREES), 
GIVEN WIND SPEED AT 1 KM (m/sec), AT KSC GEOGRAPHICAL AREA, 

FOR SPLINE CURVE INTERPOLATION (Continued) 


KSC 

m/sec 

a* 

0 

,10 

.20 

CUMULATIVE FREQUENCY 
.30 .40 .50 .60 .70 

.80 

.90 

1.00 

b** 

JUN 

0<v<5 

2 

0 

56 

90 

120 

148 

167 

197 

220 

248 

285 

360 

10 


5<v<10 

3 

0 

62 

90 

115 

141 

172 

188 

217 

241 

265 

360 

16 

II 

10<v<15 

6 

0 

93 

155 

183 

202 

216 

224 

232 

242 

257 

360 

19 


15<v 

0 

0 

104 

128 

176 

195 

210 

220 

232 

253 

280 

360 

10 

JUL 

0<v^5 

11 

0 

84 

126 

149 

172 

185 

202 

221 

241 

264 

360 

16 


5<v^l0 

17 

0 

117 

141 

160 

179 

201 

220 

237 

250 

267 

360 

16 


10<v<15 

17 

0 

107 

147 

164 

205 

220 

236 

246 

260 

274 

360 

16 


15<v 

0 

0 

170 

218 

223 

229 

232 

240 

256 

268 

290 

360 

12 

AUG 

0<v<5 

3 

0 

38 

83 

119 

138 

156 

177 

202 

240 

277 

360 

10 


5<v<10 

10 

0 

83 

115 

133 

158 

185 

210 

227 

242 

263 

360 

15 


10<v<15 

3 

0 

33 

154 

191 

210 

225 

237 

245 

254 

272 

360 

14 


15<v 

0 

0 

45 

80 

216 

221 

225 

229 

233 

242 

255 

360 

20 

SEP 

0<v<5 

. 6 

0 

45 

72 

96 

116 

137 

169 

206 

245 

301 

360 

7 


5<v<10 

5 

0 

40 

63 

81 

94 

110 

135 

168 

217 

255 

360 

16 


10<vil5 

7 

0 

44 

55 

75 

87 

99 

135 

176 

207 

251 

360 

19 : 


15<v<20 

7 

0 

26 

43 

64 

94 

146 

197 

210 

240 

262 

360 

16 


20<v<25 

3 

0 

43 

52 

59 

73 

111 

209 

215 

225 

245 

360 

16 

i 

25<v^30 

0 

0 

31 

83 

122 

138 

162 

227 

235 

260 

287 

360 

11 


30<v 

0 

0 

20 

105 

196 

206 

213 

244 

254 

290 

330 

360 

6 

OCT 

0<v<5 

2 

0 

18 

48 

73 

98 

133 

176 

223 

271 

312 

360 

6 


5<v<10 

2 

0 

20 

43 

58 

71 

90 

122 

199 

258 

301 

360 

8 


10<v<15 

1 

0 

11 

30 

55 

67 

79 

92 

114 

203 

270 

360 

14 


15< 

0 

0 

6 

16 

29 

67 

90 

219 

242 

262 

290 

360 

6 

NOV 

0<v<5 

2 

0 

18 

52 

71 

95 

151 

218 

248 

272 

295 

360 

8- 


5<v<10 

2 

0 

29 

56 

72 

89 

115 

175 

223 

256 

294 

360 

8 


10<v<15 

2 

0 

25 

47 

71 

95 

151 

218 

247 

272 

291 

360 

10 


15<v<20 

0 

0 

12 

58 

195 

227 

252 

270 

282 

287 

295 

360 

12 


20<v‘ 

0 

0 

34 

56 

214 

236 

246 

259 

268 

281 

304 

360 

8 
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Table 3-2. CONDITIONAL CUMULATIVE DISTRIBUTION OF WIND DIRECTION (DEGREES), 
GIVEN WIND SPEED AT 1 KM (m/sec), AT KSC GEOGRAPHICAL AREA, 

FOR SPLINE CURVE INTERPOLATION (Concluded) 


KSC 

m/sec 

a* 

0 

.10 

.20 

CUMULATIVE FREQUENCY 
.30 .40 .50 .60 .70 

.80 

.90 

1.00 

b** 

DEC 

0<v<5 

2 

0 

29 

94 

122 

148 

185 

218 

247 

285 

313 

360 

6 


5<v<10 

7 

0 

52 

75 

104 

142 

187 

218 

252 

281 

307 

360 

7 


10<v<15 

4 

0 

46 

84 

153 

189 

213 

240 

260 

285 

306 

360 

9 


15<v<20 

3 

0 

72 

193 

222 

239 

264 

279 

284 

290 

299 

360 

11 


20<v 

0 

0 

100 

156 

237 

253 

270 

285 

300 

315 

332 

360 

2 


* 




** 


dd) 

■ dF 
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Table 3-3. CONDITIONAL CUMULATIVE DISTRIBUTION OF WIND DIRECTION (DEGREES), GIVEN WIND 
SPEED AT 1 KM (m/sec) FOR SPLINE CURVE INTERPOLATION AT VAFB GEOGRAPHICAL 
AREA, FOR SPLINE CURVE INTERPOLATION 



— 





CUMULATIVE FREQUENCY 






VAFB 

m/sec 

a* 

0 

.10 

.20 

.30 

.40 

.50 

.60 

.70 

.80 

.90 

1.00 

b** 

JAN 

0<v^5 

0 

0 

2 

31 

86 

149 

173 

211 

251 

299 

322 

360 

6 


5<v<10 

0 

0 

0 

18 

52 

143 

180 

220 

286 

314 

325 

360 

5 


10<v^l5 

0 

0 

0 

6 

41 

143 

183 

213 

242 

308 

327 

360 

5 1 


15<v<20 

0 

0 

0 

0 

10 

137 

156 

176 

192 

202 

215 

360 

28 


20<v 

0 

0 

120 

160 

179 

191 

202 

211 

236 

260 

285 

360 

26 

FEB 

0<v<5 

0 

0 

14 

34 

97 

140 

161 

189 

237 

267 

310 

360 

7 ; 


5<v^l0 

0 

0 

0 

15 

38 

136 

162 

195 

281 

313 

335 

360 

5 


10<v<15 

0 

0 

0 

2 

24 

124 

153 

178 

232 

310 

335 

360 

4 1 


15< 

0 

0 

20 

118 

135 

150 

169 

176 

185 

198 

300 

360 

4 ' 

MAR 

0<v<5 

0 

0 

24 

82 

131 

157 

184 

216 

254 

307 

330 

360 

5 ! 


5<v<10 

0 

0 

0 

6 

22 

85 

192 

275 

308 

322 

335 

360 

4 ' 


10<v<15 

0 

0 

3 

17 

135 

164 

209 

291 

310 

330 

345 

360 

4 


15<v 

0 

0 

1 

17 

30 

138 

169 

181 

194 

208 

321 

360 

6 

APR 

0<vs5 

4 

0 

45 

79 

112 

152 

200 

227 

265 

300 

323 

360 

6 


5<v<10 

0 

0 

1 

15 

34 

159 

195 

251 

297 

315 

330 

360 

4 


10<vsl5 

0 

0 

3 

9 

21 

32 

133 

283 

305 

320 

340 

360 

4 


15<v 

0 

0 

0 

0 

0 

10 

80 

136 

170 

325 

350 

360 

4 

MAY 

0<v^5 

0 

0 

16 

52 

103 

148 

189 

223 

270 

299 

323 

360 

5 


5<v<10 

0 

0 

5 

16 

32 

101 

191 

307 

322 

330 

345 

360 

4 


10<Vil5 

0 

0 

6 

30 

305 

311 

317 

322 

328 

330 

334 

360 

4 


15<v 

0 

0 

0 

0 

0 

0 

0 

40 

305 

330 

345 

360 

2 

JUN 

0<v<5 

0 

0 

9 

48 

79 

193 

205 

266 

294 

322 

340 

360 

5 


5<v<10 

0 

0 

2 

12 

29 

86 

167 

308 

320 

330 

342 

360 

4 


10<v<15 

0 

0 

9 

20 

30 

307 

317 

326 

330 

335 

340 

360 

3 


15<v 

0 

0 

0 

0 

0 

30 

308 

314 

320 

326 

332 

360 

6 

4 

: _ d(|. 

® ■ dF . 

If . 

1 =0 



** 

b = -^ 
dF^ 
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Table 3-3. CONDITIONAL CUMULATIVE DISTRIBUTION OF WIND DIRECTION (DEGREES), GIVEN WIND 
SPEED AT 1 KM (m/sec) FOR SPLINE CURVE INTERPOLATION AT VAFB GEOGRAPHICAL 
AREA, FOR SPLINE CURVE INTERPOLATION (Concluded) 


VAFB 

m/sec 

a* 

0 

.10 

.20 

CUMULATIVE FREQUENCY 
.30 .40 .50 .60 .70 

.80 

.90 

1.00 

5 ** 

JUL 

0<v<5 

0 

0 

10 

38 

76 

140 

176 

213 

259 

305 

331 

360 

5 

5<v^l0 

0 

0 

6 

16 

30 

56 

157 

295 

313 

324 

332 

360 

5 


10<v 

0 

0 

0 

10 

310 

315 

320 

325 

330 

335 

340 

360 

4 

AUG 

0<v<5 

0 

0 

15 

47 

115 

164 

196 

255 

284 

306 

330 

360 

5 


5<v^lO 

0 

0 

3 

19 

44 

142 

190 

308 

319 

328 

340 

360 

3 


10<v<15 

0 

0 

0 

0 

131 

149 

170 

315 

325 

335 

345 

360 

1 


15<v 

0 

0 

0 

0 

131 

149 

170 

315 

325 

335 

345 

360 

1 

SEP 

0<v<5 

0 

0 

11 

53 

101 

133 

159 

181 

216 

280 

322 

360 

5 


5<v<10 

0 

0 

6 

26 

60 

124 

146 

161 

281 

330 

345 

360 

3 


10<v" 

0 

0 

0 

0 

124 

133 

140 

330 

335 

340 

345 

360 

1 

OCT 

0<v<5 

2 

0 

41 

86 

124 

145 

164 

191 

232 

268 

311 

360 

6 


5<v<10 

0 

0 

0 

6 

14 

27 

47 

148 

200 

287 

330 

360 

3 


10<v 

0 

0 

0 

3 

7 

13 

27 

32 

330 

340 

350 

360 

1 

NOV 

0<v^5 

0 

0 

12 

43 

115 

149 

168 

206 

256 

295 

325 

360 

4 

5<v^l0 

0 

0 

5 

28 

115 

136 

155 

179 

217 

266 

316 

360 

6 


10<v^l5 

0 

0 

0 

1 

11 

25 

131 

153 

182 

210 

240 

360 

8 


15<v 

0 

0 

9 

31 

142 

161 

172 

178 

185 

194 

270 

360 

2 

DEC 

0<v<5 

7 

0 

61 

97 

140 

162 

183 

214 

269 

293 

320 

360 

5 


5<v<10 

0 

0 

0 

9 

30 

124 

162 

205 

262 

306 

330 

360 

2 


10<v<15 

0 

0 

0 

13 

58 

146 

212 

277 

312 

328 

340 

360 

3 


15<v 

0 

0 

11 

169 

180 

191 

221 

228 

236 

260 

326 

360 

3 


* 
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Section IV 


LONG-CRESTED OCEAN WAVE MODEL 


4.1 GENERAL 


Long-crested wave heights n can be represented as a two-dimensional 
Fourier integral. 


n(t,x) 





fi(a,k) 


gi (kx-at) 


dk da 


(4.1.1) 


where n is the Fourier spectrum, t is the time coordinate, and x is the distance 
coordinate taken in the direction of travel of the long-crested waves. The term 
k is the radian wave number and a is the radian frequency. The inverse of 
Equation (4.1.1) is given by 


n(a,k) 


+0O 

j n(t,x) 

~00 


-i (kx-at) 
e 


dx dt 


(4.1.2) 


For deep water and shallow waves, the dispersion relationship (see Reference 7) 


2 

a 




(4.1.3) 


is commonly used, where g is the acceleration of gravity. This relationship 
relates the wave number to the frequency so that waves with small wave lengths 
will also have small cycle times. Since a and k are directly related, then 

yN ^2 

n(a,k)=2TT r)(cr) 6(a -k) , where 6 is the Dirac delta function. Then Equation 
4.1.1 can be written as 


n(t,x) 


2-n 



n(a) e 


i (kx-at) 


da 


— 00 


(4.1.4) 


where k is now a function of a as given by Equation (4.1.3). The above expres- 
sion can also be written as 


p-l-oo 

n^(a,x) e ^^^da (4.1.5) 

— 00 

4-1 


n(t,x) = ^ 


where 


= ti(ct) e 


ikx 


(4.1.6) 


The inverse Fourier transform is then given by 

+00 


n^(a,x) = 


n(t,x) dt 


(4.1.7) 


Notice that it would have been possible to retain k-dependence instead of the 
a-dependence. Then in Equation (4.1.4) the result would have been 


+00 


n(t,x) = 


27T 


" V i(kx - at) 
n(k) e ^ dk 


(4.1.8) 


which could have been written as 

+00 


n(t,x) = 


ikx ,, 
e dk 


(4.1.9) 


Thus the present calculations can be carried out either in frequency or 
wave number space. 


4.2 RELATIONSHIP OF FOURIER SPECTRUM TO POWER SPECTRUM 

To relate the Fourier spectrum to the power spectrum, following Reference 5, 
for a specific value of x the following equation can be written 


T 


R (t,x) 

nn 


lim 1 
T-^ T 


n(t,x) n(t+T,x) dt 


(4.2.1) 


T 

2 


The Fourier transform (FT) of the above expression is 

-Ho 




FT[R (t,x)] = ^ I 

nn ’ ' ^ T-^ T J 


n(t,x) 


T 

2 


n(t-fx,x) 6^°^^ dt 


dt (4.2.2) 


4-2 




FT[R^^(t,x)] = I n(t,x) rijj ( 0 ,x) dt 

T 

“ 1 


lira 1 '•* 
T->oo T *^x 


(o,x) 




(4,2,2) 


Taking ensemble averages and using Equation (4.1,6), the power spectrum 
S(a) is given by 

S(0) = FT[R^^] = FT[<R^^>] = (4.2.3) 

where R is the autocorrelation, 

nn 

R^^(t) = <n(t) n(t + t)> (4.2.4) 


4.3 OCEAN WAVE SPECTRA FOR LONG-CRESTED WAVES 

One ocean wave spectral distribution which is commonly used is the Pierson- 
Moskowitz model (Reference 6). This model applies to a ”fully aroused sea." 

This is a sea in which all the Fourier components have been excited to the fully 
developed state by the wind stress. To attain these conditions, a wind must 
blow for a sufficiently long time over a sufficiently long fetch (distance over 
unobstructed water). The statistical characteristics of the "fully aroused sea" 
are considered constant. 


The Pierson-Moskowitz spectrum can be written as 


2 

S(a) = 7T a exp 
a 



where a and 3 are nondimensional universal constants 
a = 8.10 X 10“^ 

3 = 0.74 

and g is the acceleration of gravity given by 

-2 

g = 9.8 (m - sec ) 


(4.3,1) 


4-3 



V is the mean wind at the 19.5 m reference level in m-sec , and a is the 

o , 2 -1 

radian frequency in radian-sec . The units of the spectrum S are m /rad sec 
Since the autocorrelation is the inverse Fourier transfoinn, (FT , of the power 


spectrum^ 


R (o) = FT-^S] =1^ ^LO|_exp -a 


(4.3.2) 


4.4 NORMALIZATION OF THE VARIABLES 

To simplify the power spectrum of Equation (4.3.1) the frequency is 


normalized to be 


(4.4.1) 


Since at must be dimensionless, the dimensionless time must also be defined as 


T =t('^ 


(4.4.2) 


so that 


o)T = at . 


(4.4.3) 


For this case the spectrum becomes 


S(a)da = S(a) 3 ^ do) = S(a))d 60 . 
dco 


(4.4.4) 


Equation (4.4.4) can be written as 

Trav ^ -(o))“^ 

o e 

S(o))da) = — ^ do) 
g 6 03 


47T a^ do3 

n 5 

CO 


(4.4.5) 


From Equation (4.1.3) 


k = 2 .^= 

g „ 2 


C4.4.6) 


4-4 



Thus, k can be normalized to be 


2 

K = = 0) 


1/2 


gg 


Since kx must be dimensionless, 




1/2 


X = X 


V 


2 


(4.4.7) 


(4.4,8) 


Then normalizing. 


P = 


(4.4.9) 


Then from Equation (4.1.4) 

+0O 


p(t,x) = ^ 


P (w) e 


i(KX ” OOT 


+ 0O 


e dw . (4.4.10) 


Then the autocorrelation and spectrum become 


R 


nn 


2 PP 


= R = <p(t) p (t + At)> 


(4.4.11) 


and the power spectrum becomes 


5 ( 03 ) ^ 
a 


l^ = ,i;(oo)=^ f <p*(o)) p(io)> 

0) 




lim 1 


(4.4.12) 


The power spectrum ^|^(o3) is plotted on Figure 4-1. It can be seen that it is a 
narrow band spectrum with a maximum at 0)=0.946. 
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I 



5 A 



Figure 4-1. DIMENSIONLESS POWER SPECTRUM ip 
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Section V 


MONTE CARLO SIMULATION OF LONG-CRESTED WAVES 


5.1 GENERAL 

Two possible methods of simulating ocean waves will be presented in this 
section. One method uses a discrete Fourier series with randomly chosen 
coefficients and phase angles. The second method is a control system approach 
in which white noise is inputted through a linear system, and the output signal 
has the same statistical behavior as the simulated wave motion. 

5.2 DISCRETE FOURIER SERIES SIMULATION 

The dimensionless wave height over a period of dimensionless time 
for a given value of x can be represented by a Fourier series (Reference 5) 


P (t,x) 


N 

2 

= I An(u,x) 

N 


e 


(-io) t) 
n 


(5.2.1) 


where 


and 


0 ) = = nAoj 

n T 

P 



(5.2.2) 


and is given by 


T 


A = — 
n T 


j lim 1 f V 

p e ^ dr = ^ — (p (oj ) 

^ T ^ T X n 

P P 


P (tOn) e i 


T 

_ ^ 
2 


(5.2.3) 


then 


. (i<„x) Ixm 1 , ss 

A=ae a= ^“(pC^jJ)) 

n n nT^T n 

P P 


(5.2.4) 


5-1 


I 


Then 


I 


4^ 

„ _ V „ X-w t) (5.2.5) 

P Z, a e n n 

1 .- N 

n- 2 


This can be expanded to 




p = ^ (a + ia „)(cos(< x-w t)+1 sin(K y-“ t)) 

^£V)il XjXl ilLl 1111 

" ~1 


(5.2.6) 


Since p is real, the imaginary term must be zero. As shown in Reference 5, 

this is because is an even function and is an odd function about 

R,n I,n 

03^ = 0. Then by expanding Equation (5.2.6) 


^ 2 

p = y cos(k X“^ 

“ R,n n"^ n 


sin(K Y-o) t) 
I,n n"^ n 


(5.2.7) 


Since is even and cos(k v-oj t) is even, the product is even. Similar! 

R,n ^ n^ n ' ^ 

since a, and sin(K y-co t) are both odd functions then this product is an even 
I,n n n 

function. Thus, the above can be written as 


N/2 

P = “r.o 2 I [«R,n cos(k^x-V> - “l.n sln(K^x-v) • -2 .8) 

n=l 


It has been shown in Reference 5 that 
where 


and are random Gaussian variable 

R,n I,n 




1 

2 

27T<a > 
n 



+ a 


2 

I,n 


)/2<a;> 


(5.2.9) 


5-2 



and 


<a„ > = 

R,n 


<a-r „> = 0 

I,n 


(5.2.10) 


Furthermore, 





> 


From Equation (5.2.4), the following is obtained 


(5.2.11) 


<a a > = 
n n 


2 

<a^ > 

R,n 


<a^ > 

I,n 


= 2<a > = 
n 


lim 

X ->oo 

p 


<p (o) ) p*(0J )> 
n n 


(5.2.12) 


Then, from Equation (4.4.12) 


2 

2<a > 
n 


lim 

X 

P 


T 

P 


lim 

X 

P 


ip(u)n) when n ^ 0 


(5.2.13) 


However, when oj =0, 
n 


<a^ >=0, the above expression becomes 

I,n ^ 


2 

<a > 
o 


lim i^(o) lim , , . Aoj 
— -- - lijCo) XT” 

X-X30 X T ^ 2 n 

p p p 


where ipCix)^) is given by Equation (4.4.12). 


(5.2.14) 


If the transformations 


= ot cos e 
R,n p,n n 


ot 


I,n 


a sxn £ 
p ,n 


are used, then Equation (5.2.8) becomes 


(5.2.15) 


N/2 

p = a cos £ + 2 y a cos(k X"to t+£ ) (5.2.16) 

p,o o ^-P,n nnn 

n=l 


5-3 


where the joint distribution of a and e are obtained by transforming 

np n 

Equation (5.2.9) to obtain 


f (a 


which becomes 


(5.2.17) 


f(a )da 
p,n p,n 



-ot 


p,n 


/ 2<ot.^> 


h da 


p,n 


(5.2.18) 


and 


- If ■“'n • 

Equations for randomly sampled values from Equation (5.2.18) and Equation 
(5.2.19), as shown in Reference 5, can be derived by using the uniform random 
number R in the following equations 


e 


n 


2ir R 

en 


(5.2.20) 


a = [2<aJ> In (^) ] 

" \,n 


1/2 


(5.2.21) 


2 

where is given by Equations (5.2.13) and (5.2.14). 

Equation (5.2.16) has been given in a number of References (for example. 

Reference 9) as a probablistic model for ocean waves^ but it was not derived 

from basic equations as in the present case nor was it used as a simulation 

procedure. It had also been often assumed that a was constant rather than 

p ,n 

a random variable. 
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5.3 CONTROL SYSTEM SIMULATION 

In the control system simulation, a random white noise signal, I, is 
inputted into a control system, H, designed so that the output p(x»t) has the 
desired statistical behavior (Figure 5-1). Then using Equation (4.4.10) 


(«< 3 »X) = p(u)e^'^^ = H(o)) 


(5.3.1) 


where 


FT ^[p^(o),x)] = pCtjX) 


It can be seen that 


(5.3.2) 


“^Py p > = <p(0)) p (w)> = H(u) H (to) <I(co) I ((0)> . 


X ' X 

This can be written as 


(5.3.3) 


ip(iti) = H(to) H (to) )|)T.(a)) 


(5.3.4) 


Since the input signal is white Gaussian noise with a unit spectrum, the wave 
spectrum is given by 


<p(co) = H((o) H (to) 


(5.3.5) 


In normalized form it becomes 




47re 


- 0 ) 


-4 


= H(o)) H (o)) 


(5.3.6) 


The system function H(a>) which will transform the noise signal can be written 
in terms of a Fourier spectrum A(o)) and a phase angle 4)(o)) 


H((o) = A(to) 


(5.3.7) 
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Letting cf)(u)) be zero, obtain 


H(u) = A(to) = . 

Then the output signal can be expressed as 

P = I(o)) . 

X 

If the impulse response function h is 

h(r;x) = Ft"^ ; 

then the wave signal is expressed as 

p(t;x) = (p ) = h(x;x) * i(t) 

X 

where * refers to the convolution integral 


p(t;x) 


h(T';x) I(t-t') dt’ 

— 00 


(5.3.8) 


(5.3.9) 


(5.3.10) 


(5.3.11) 


(5.3.12) 


Thus, by generating a noise signal and convolving it with the appropriate 
system function, a wave field in time and position can be generated. If a 
signal is desired for a point moving along the surface at some velocity v, 
then x~vt; and then equation (5.3.12) will give the wave signal in time for 
an observer moving with velocity v. 
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Section VI 

WAVE CHARACTERISTICS 


6.1 GENERAL 

The quantities that describe the wave characteristics, specifically the 
probability distribution of wave heights, wave slopes, and velocities, are 
given in this section. 


6.2 PROBABILITY DISTRIBUTION OF WAVE HEIGHTS 

As discussed in References 7 and 10, the wave height has a Gaussian 
probability distribution given by 


f(p) = 


T/2 


-2 — ^ = \rro 


( 6 . 2 , 1 ) 


2 

where the N(p;<p> <p >) refers to a normal distribution of p with a mean of 

2 2 

<p> equal 0 and variance <p >. The <p > is given by 


+00 


<P^> 


2tt j 


i|j(a))du) = M 


0,0 




-4 


= 4 ^ dco = 1 . 


( 6 . 2 . 2 ) 


Thus, the normalized wave height has a Gaussian distribution with a mean of o 
and a variance of 1. The M. . are moments of the spectrum and will be used more 
extensively later. 


« 1 1 

«i.j - 2? J “ 


k((o)^ y(o)) do) 


Since very small waves will have very little effect on the Space Shuttle 
Solid Rocket Booster, it is worthwhile to truncate the spectrum, at some maximum 
frequency, This will not change the mean value of the wave height distri- 

bution which will remain 0; however, the variance of the wave height will be 
T 

reduced to ^ which is the integral of the power spectrum up to 

“o.o = ¥ j I|^(a3)dw = e (6.2.3) 

o 
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The distribution of wave height will now be given for the truncated spectrum 
case by 


f^(p) = N(p;o,M^^^) 


(6.2.4) 


The values of M for different values of oi are given in Table 6*~1, 
o, o max ® 


Table 6-1. VALUES OF MOMENTS OF POWER SPECTRUM TRUNCATED AT VARIOUS VALUES 
OF to 

max 


03 

max 

0,0 

0,2 

2,0 

”1.1 

00 

1.0 

00 

1.772 

-3.625 

10 

0.9999 

8.63 

0.876 

-3.226 

3.16 

0.99 

4.04 

0.787 

-2.362 

2.10 

0.95 

3.49 

0.669 

-1.740 

1.75 

0.90 

1.775 

0.572 

-1.387 


It is of interest to check the distribution of wave heights given by the Fouriei 
series simulation procedure, Equation (5.2.11). The central limit theorem 
states that a random variable made up of a random sum will have a Gaussian 
distribution, permitting 

f(p) = N(p;<p>,<p^>) . (6.2.5) 

It can be seen that by taking the ensemble average of Equation (5.2.16), where 
use is made of the fact that a and e are statistically independent^ 

N/2 

<p> = <a xcose > -f 2 7 <a > <cos(k X"^ 9+^ )> - (6.2.6) 

p,o o ^ P,n ^ n^ n n^ 

n=l 
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Using Equation (5.2.19), obtain 

2tt 

I 

<cos e 


= [ 


de 


<cos(k )> 

^ n^ n n' 


o 

r 2 tT 


COS TT— = 0 
o 2 tt 


de 

cos(<nX-%®+^n) IT 


= 0 


(6.2.7) 


( 6 . 2 . 8 ) 


so that <p> = 0. For <p >, there is 


2 2 2 N/2 2 2 

<p > = <a^ ^><cos e^> + 4 ^ <a^ ^><cos ^ ' (6.2.9) 

* n=l * 

Cross-product terms are neglected since a and e are statistically independent. 
Using Equation (5.2.18)^ obtain 


2 

<a > 
p,n 


a 


2 Ipn 


n 


e““ p,n/2<a^>) da = 2<a^> 
n / p ,n n 


(6.2.10) 


Similarly, 


2 1 
<cos (k 0+^ )>= 

n n 2 tt 


2tt 


2 1 
cos (k x-w 0+e^) de = ^ 
n n n n 2 


( 6 . 2 . 11 ) 


Therefore. 


2 2 2 1 
<p > = <a > + 4 y <a > = -TT- [i|^(o)Ato+2 T )Ao)] 
o n Ztt n 

n=l n=l 


( 6 . 2 . 12 ) 


OO 


27T 


J 


— CO 


\|;(ai)dw = M 

0,0 


This gives the same result obtained in Equation (6.2.1). It is also of interest 
to check the control system simulation to see if it gives the appropriate distri- 
bution of wave heights. From Equation (5.3.1) 
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P(T,X) 


(6.2.13) 


+00 

I 


H(w) I(o)) e do) 


Inputting a white noise Gaussian process into a linear system gives 
output so that 

2 

f(p) = N(p;<p>,<p >) 


where 


<p> = FT ^ H(to) <I(o))>] = 0 


o 2 1 

and <P > = ■^ 


r +~ 


i|/(to)do) = 


4od 


H(o)) H*(a))do) = M 


0,0 > 


the result agrees with the previous result. 

6.3 PROBABILITY DISTRIBUTION OF WAVE SLOPES 

The normalized wave slope is given by 


dp _ 1 

dx 2 tt J 


^oo 


p (o)) (iK) e^^^^ do) 


As pointed out previously, a Gaussian random signal passing through 
system will remain Gaussian; therefore ^ 

f(p^) = N(p^;<p^>,<pJ>) . 

From Equation (6.2.13) it can be seen that 

(Ik) H(to) I(t»)) e ^^^dco 



Gaussian 

(6.2.14) 

(6.2.15) 

(6.2.16) 


(6.3.1) 
linear 

(6.3.2) 

(6.3.3) 
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Then <p > = 0 as before, and the spectrum is given by 
X 

<p p > = H(co) H*(m) = tc^ijj(a}) 

X X 


(6.3.4) 


2 1 2 

<p > = ~ K 

X 2tt J 


K i(;(a))d{jo == 2 


(6,3.5) 


where, using the dispersion relationship (Equation (4. 4 .7)) 


( 4 e ^ 

M ^ = 4 oj — do3 . 


Letting y = — the above becomes 


CO —y 

- dy . 

y 


As before, interest lies in the slopes due to the larger waves^ and the advan- 
tage of removing the effects of the smaller waves which will not affect the 
Space Shuttle Solid Rocket Booster. By retaining the lower frequency fraction 

of the spectrum and letting the spectrum be zero for values greater than co » 

max^ 

the slope distribution is now given by 




where 


0,2 TT 


i|<(a))daj =1 - dy = El (-| — ) 

^ max 

max 


and El refers to the exponential Integral. This is evaluated for several 
truncations, and the values are shown in Table 6-1. 
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6.4 FREQUENCY DISTRIBUTION OF WAVE VERTICAL VELOCITIES 

The normalized vertical wave velocity is given by 


+00 

dp 1 f ^/ \ / . N i(icx~“T) , , 

^ = Pt “ ^ J p (W) (- 103 ) e dm 

— 00 

As previously discussed, will have a Gaussian distribution 


f(p^) = N(p^;0 ,M2^^) 

where the variance is given by 


(6.4.1) 


(6.4.2) 


M, 


2,0 


2 

<p^> 


= - f 

2ir J 


0 ) if/(m)do) = 4 


OO —4 

f 2 e”“ 
m — j dm 

J /i\ 


= JiT^ 1.772 


(6.4.3) 


Primary interest is in the vertical velocities of the large waves. If the 
higher frequencies are removed from the spectrum^ the result will be a truncated 
wave vertical velocity distribution 




(6.4.4) 


where 



1 


7T 


0 ) 

• max 2 

0) ip(a))do3 = ^ 

o 


0 ) -"+ 
max 

\ dm 

J m 


-y^ . ^ .1 . 

2 dy =''-^ erfc (- 5 -) 


max 


2 

03 

max 

and erfc refers to the complementary error function integral, 

are given in Table 6-1 for various values of o) 

max 


(6.4.5) 


These results 
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6.5 FREQUENCY DISTRIBUTION OF HORIZONTAL WAVE VELOCITY 

Consider a point on the wave surface at some height n * time t, and position 
X. At some later incremental time the wave surface will move to a new point 
X+^X 3.t time t+At but at the same height n. 


dp = + (|7)dT = 0 . 


The normalized horizontal velocity is given by 


(6.5.1) 




dx 


9t 

IP 

9X> 


(6.5.2) 


The joint Gaussian distribution of and p can be expressed as 


Px “ ^ ^1.1 ^X^ “o,2 


where 


" = ^2,0 V2 - ^ 1,1 • 


This joint Gaussian can be transformed to the new variable 


(6.5.3) 




to give 


f(x_.y) = ~^i/2 |d| exp [-y^ (-^ + — ^ + M 

\ xj X^ ’7 


C6.5.4) 


where the Jacobian is given by 

iJi =1^ 
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This can be integrated over the y variable from -«> to -H» to give 


> 1/2 


f (X,) = -f 


M o + 2 

0,2 


^ 1,1 ^ ^ ^ 2,0 


(6.5.5) 


This gives the distribution of the horizontal velocity. 


It can be found that the maximum or mode of the distribution is given by 


M 

T 


M, 


1,1 


M 


0,2 


(6.5.6) 


Then 


f(x^) 


> 1/2 


■n’M _ Cx^ “ 2x X^ + w ^ ° ] 
o2 M_ 


0,2 


> 1/2 


«o,2 X,) * ^2 J 

0,2 0 , 2 -^ 


(6.5.7) 


which can be expressed as 


> 1/2 


f (Xt-) 


),2 ^ 


ttM I (X_-X^)^ + ~f 

0.21 T T 

o 


, 2 ^ 


(6.5.8) 


M 

This is a Cauchy distribution with a maximum at x = i the 

M 

distribution is symmetrical at • 
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For purely progressive waves, since k and o) must be of opposite signs, and 

2 

remembering that k - o) , 



CO 

Kajip(a))d(jj = -4 
o 



3 e 


- 0 ) 


-4 


0) 


do) 


OJ 


(6.5.9) 


This can be rewritten with y = — ^ as 


0) 


M. 


1,1 


■^4 dy = - r(f) = -3.625 


(6.5.10) 


where T represents the gamma integral. 

As before, the Interest is in the horizontal velocity of the larger waves. 

The spectrum will extend only to and the integral becomes the incomplete 

Gamma function 




03 

max 




max 


(6,5.11) 


6.6 VELOCITY FIELD BELOW THE WATER SURFACE 

As discussed in Reference 9, the horizontal and vertical velocity below 
the surface can be related to the motion on the surface by using the referenced 
'*Long-Grested Gaussian Linear Eularian Model of Random Waves." If Equation 
5.2.16 is used for the surface wave model, then the horizontal velocity u and 
the vertical velocity w at height z is given by 


and 


N/2 

u = a cos £ + 2 J a (exp(-K z)) cos(k 0+e ) (6.5.12) 

p,o o p^n n n n n 


w 


N/2 

a cos e + 2 7 a (exp(-K z)) sin(tc y~ 03 0+e ) (6.5.13) 

p,o o p,n n n'^ n n^ ^ 


where Z gives the depth below the surface. 
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Section Vi I 


CONCLUSIONS 

This report presents the pertinent features of a study of the effect on 
the Space Shuttle Solid Rocket Booster of oceanic impact and the subsequent 
recovery. Data on the statistical distribution of ocean currents in the 
planned splashdown areas are given in a form suitable for use in simulation 
studies. This is because spline fitting can be used for interpolation in the 
random sampling procedure of the simulation. 

The influence of 1-kilometer winds on wave behavior is also discussed. 
The ocean model based on random function analysis and assumed long-crested 
waves is developed. Using the Pierson-Moskowitz power spectrum distribution 
of wave height, wave slope and wave velocity are given. The velocity of the 
fluid below the surface based on the linear Eularian model is also shown. 

Two simulation procedures are developed, but no extensive numerical 
calculations are presented. Such numerical calculations for the simulation 
models as well as extension of the results of the long-crested wave model 
are currently being done for the case where waves can move in all directions. 
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